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Abstract. - The observed discontinuity in strain-rate for a two-dimensional foam undergoing 
shear in a (circular) Couette system is explained in terms of the continuum (Herschel-Bulkley) 
model. It is attributed to the finite difference between yield and limit stress. 
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j Introduction.— Many of the existing observations in 
2D foam rheology [1-4] may be accounted for within the 
general continuum model that incorporates bulk dissipa- 
tion and, where appropriate, wall drag [5-8]. One feature 
that has so far escaped such an explanation is the discon- 
tinuity of the derivative of velocity (or of the shear rate), 
reported by the group of M. Dennin [3,9] for the case of 
a 2D foam in a circular (Couette) rheometer. Figure 1 
shows the data. These measurements were made using 
the variety of 2D foam that is also known as the Bragg 
raft. 

This derivative in question is everywhere continuous in 
the existing continuum model. In the present paper we 
show that a natural refinement of that model gives rise 
to just what is observed, and our simulations of the effect 
offer opportunities for further experimental tests. 

The refinement that is required is the incorporation of a 
stress-strain relation that has a peak or overshoot before 
it settles down to its asymptotic value. That is, there is a 
yield stress that is in excess of the eventual limit stress, and 
it is the difference between these quantities that relates to 
the discontinuity (see Fig. 3). This has for long been a 
familiar feature of measurements and simulations of foam 
rheology (e.g. [11]). 

However, as elsewhere in the physics of foams, the out- 
come depends on the detailed history (and even prehis- 
tory) of the experiment, so it is essential to study simula- 
tions of the velocity profile in various cases. We shall see 
that this discontinuity does indeed depend on the experi- 
mental protocol. 
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Fig. 1: Reprinted figure 5 from [10]. In the case of circular 
(Couette) geometry, the scaled azimuthal steady-state velocity 
v(r)/(rfl), where r is the radial position and Q is the rotation 
rate of the outer cylinder, may become constant at some in- 
ternal point (see also figure 2). In the data of [3,9, 10], taken 
for a moving outer boundary and a fixed inner one, there is a 
discontinuity of d(v(r)/(rQ))/dr at this point. Note that this 
published figure uses the symbol v(r) to denote v(r)/(rfl) in 
the notation of the present paper. The solid and dashed lines 
show least-square fits to power-law and exponential, respec- 
tively. This figure is used with permission from M. Dennin [10], 
Copyright 2008 by IOP Publishing Ltd. 
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Fig. 2: Sketch of velocity profile for a moving outer boundary 
in the case where continuous shear is confined to R < r < r c . 



The effect should not, of itself, be restricted to the 2D 
foam system in question, but we are not aware of any 
reference to it in 3D foams. In any case, it offers further 
support for the continuum model in the present context. 

In most of what we present, the original primitive ver- 
sion of the continuum model is used [5]. In this, linear 
forms are used for bulk dissipation and wall drag terms. 
Later we will indicate the straightforward modification of 
the main result that arises from the introduction of non- 
linear forms, corresponding to the later and more realistic 
development of the model [8] . 

The continuum model with distinct yield and 
limit stresses. — Clancy et al. [6] have provided an 
extensive analysis of the 2D foam continuum model for 
circular (Couette) geometry. It is a much richer scenario 
than that of simple shear between straight boundaries, 
although the two are closely similar in some limits. Local- 
isation due to wall drag is combined with geometric effects 
in this case, so that there are several qualitatively distinct 
forms of the velocity profile, according to the dimensions 
of the system, the parameters of the model, the imposed 
velocity, and whether it is applied to the inner or outer 
boundary. 

A particular feature of the circular case is the possibility 
of the existence of a range of radial positions , r, at which 
the strain rate is zero (implying rigid body rotation), as 
in Figure 1 and the sketch in Figure 2. 

The present topic relates to just this feature, specifically 
to the boundary point r = r c between the ranges of zero 
and non-zero shear rate. 

We follow closely the formalism and algorithms of the 
previous work [5-8] . In the analysis presented here we will 
restrict ourselves to a moving outer cylinder to match the 
experimental conditions of Figure 1. The local strain rate 
y(r) will thus be positive throughout. 

We begin with an elementary steady-state treatment 
that yields a formula for the discontinuity, before embark- 
ing on a time-dependent simulation. 

The constitutive relation expressing local stress <r(r) as 



Fig. 3: Plot of the stress-strain function f(x; y, z) of eqn. 2 
where we have set y = — = 1.08 and z = — = 2 for illustra- 
tive purposes. Numerical evaluation of eqn. 3 then results in 
(3 = 4.1076. The form of f(x;y,z) as given by eqn. 2 is only 
valid for the case where the strain rate is always positive. 



a function of radial position r, local strain 7(7") and local 
strain rate j(r) is given by [5-8] 



VLf{l/lY\<JY/<JL,lL/-fY) + ?/7 



(1) 



where we call the first term on the right-hand side the 
elastic-plastic stress. Yield and limit stress are denoted 
by try and ox, respectively, 7^ and 7y are yield and limit 
strain. The viscosity component of stress (also called con- 
sistency) is denoted by rj. As a function / representing 
the scaled quasi-static stress-strain relation we have arbi- 
trarily chosen 

x for x < 1 

1 + 0-!) (jE^Y forl<x<z (2) 
1 for x > z 

where (3 is fixed in terms of y and z by the value of the 
maximum of / (w.r.t. x) for 1 < x < z, 



= V —l 

(1 + /?) 1 +' 3 z-1 



(3) 



See Fig. 3 for a plot of the yielding curve used in our 
simulations in our time-dependent simulations, described 
below. 

As the strain-rate discontinuity is reported for bubble 
rafts [3,9], the expression for the wall drag in the contin- 
uum model [5] is omitted. However, the key result of this 
paper would not be affected by it. 

Following from appendix A. 3 of [6] we thus have (in the 
absence of inertia) 



r 2 dr 







(4) 
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and so the radial variation of stress is given by 

A 



(5) 



where A is a constant, fixed, for example, by the stress on 
the moving boundary. 

The matching conditions for the steady state profile at 
the boundary point r c are as follows. 

• Velocity v(r) is continuous '(otherwise the strain rate 
would be infinite). 

• Total stress a is continuous since it is given by eqn. 5. 

The second condition is most relevant here. On the side 
of r c that is undergoing continuous shear with rate 7(f), 
Equation (1) reduces to: 



a = cfl + rff, 



(6) 



where <tl is the limit stress (Fig. 3). 

On the other side of r c , there is no continuous shear, so 
the total stress is given by the elastic-plastic stress (first 
term in Equation 1). Prima facie, it appears that the value 
of the elastic-plastic stress at r c can be anywhere in the 
range < a < ay, and we here assume the value a = ay 
is appropriate. The next section describes time-dependent 
simulations that validate this choice for the usual experi- 
mental procedure. 

Equating the total stress on the two sides of r c , we ob- 
tain 



7(r c ) = (oy - (T L )/rj. 



(7) 



Since the yield stress exceeds the limit stress, ay > <jl, 
the strain-rate on the shearing side is finite at r c , leading 
to a discontinuity in the strain-rate at that point, hence to 
a discontinuity in the azimuthal velocity gradient at this 
point. 

Range of possible magnitudes of the disconti- 
nuity. — The argument above allows any value of the 
elastic-plastic stress at r c that is less than ay, however, 
the possibility of having a < &l may be dismissed as fol- 
lows. 

Such a value would imply, through Equation 6, a nega- 
tive value of 7 for values of radii less than, but close to r c . 
Assuming continuity of the velocity, this in term implies 
a velocity maximum in this range, at some value r rnax . 
There is no force on an element at that radius due to the 
elastic-plastic stress (in the shearing region the elastic- 
plastic stress must be the limit stress everywhere, because 
we are dealing with the steady state), but there is a force 

present due to 8 g^i ■ Since this cannot be balanced, a 
value of a which is less than uj, is thus not possible. 




Fig. 4: Exemplary calculations for the continuum model in the 
circular geometry of figure 2, with L — 1.257? and the stress- 
strain function of figure 3. (a) shows the steady state for both 
local v(r)/r and strain-rate 7 (both expressed in dimensionless 
units) for an outer cylinder rotating at a rate fii = ^jr^ ~ 
0.1. The position r c /R ~ 1.19, where the strain-rate has a 
discontinuity separates the shearing (r < r c ) from the non- 
shearing region (r > r c ). A decrease of the rotation rate of 
the outer cylinder to Q2 = 0.08 eventually leads to the steady- 
state profiles shown in (b). Note that the discontinuity in the 
strain-rate has now vanished. 



Time-dependent simulations. To explore this 
further, we resort to the full implementation of the time- 
dependent model [6] , and simulate an experiment in which 
the system starts from static stress-free equilibrium, and 
the outer boundary is set in motion with a fixed velocity. 

Figure 4 (a) shows angular velocity and strain-rate once 
a steady state is reached 1 . The strain-rate is discontin- 
uous at r c . Also, in agreement with our argument above 
and Eqn. 7, it remains positive in the vicinity of r c , and is 
given by j(r c ) / '{(Jlv) = 0-08 for the choice of try = I.O80X 
in the simulations shown. Upon decrease of the velocity 
of the outer boundary a new steady state is reached with 
velocity and strain-rate profiles as shown in Figure 4 (b). 
Our simulations show the absence of the discontinuity of 
the strain rate for this scenario. 

As in previous work [6] our algorithm for the imple- 
mentation of the continuum model is not stable unless al- 
lowance is made for the possibility of negative strain rate 
7, even though this would appear not to occur in the solu- 
tion represented by the numerical results. We still do not 
have any explanation of this technicality. 

General continuum model. In the generalisation 
of the continuum model we introduce a non-linear depen- 
dency of the strain rate dependent term in the constitutive 



1 Movies of the simulations can be downloaded at 
http: / / www. ted. ie / physics / foams / publications . php 
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relation of Equation (1), 

<y = VLf{i/iY]VY /ol,il/iy) + v'i a , (8) 

where the exponent a is called the Herschel-Bulkcly ex- 
ponent [8]. Using the same argument as for the linear 
case, we again find a discontinuity in the strain rate at the 
boundary point r c between the ranges of zero and non-zero 
strain-rate. The value of strain rate at this point is given 
by 



7(r c ) = 



cry — o~L 
V' 



l/a 



(9) 



Conclusion. — The overshoot that plays a central role 
in the above arguments is quite familiar in foam physics, 
but is nevertheless not understood in any detailed sense. 
It is attributable to the change of structure under imposed 
continuous shear [12]. 

The meaning of this is not entirely self-evident. Suppose 
we shear a foam considerably, and then reduce stress to 
zero. Of course, the new structure is microscopically dif- 
ferent from the original one, but it is also different in an av- 
erage, statistical sense. It will, for example, be anisotropic, 
even if the original sample was not. 

The explanation of the shear rate discontinuity, as pre- 
sented here, may not be unique but seems compelling. It 
invites further tests, by varying the experimental param- 
eters and, in particular, following the protocol which is 
predicted to eliminate the discontinuity. 

This particular effect is not confined to two-dimensional 
samples. Accordingly, it may well have been noted in 3d 
rheometry at some earlier stages, but we are unaware of 
any direct such observations. In the present context, its 
explanation lends further support to the continuum model 
as the natural "mean field" model to be applied to shear 
localisation in 2D foams, at least as a first approximation. 
It also offers a test for further experiments, to verify that 
the discontinuity may be eliminated by a different proto- 
col, as seen in our time-dependent simulations. 
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